Probing a nonequilibrium Einstein Relation in an aging colloidal glass 
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We present a direct experimental measurement of an effective temperature in a colloidal glass of 
Laponite, using a micrometric bead as a thermometer. The nonequilibrium fluctuation-dissipation 
relation, in the particular form of a modified Einstein relation, is investigated with diffusion and 
mobility measurements of the bead embedded in the glass. We observe an unusual non-monotonic 
behavior of the effective temperature : starting from the bath temperature, it is found to increase 
up to a maximum value, and then decreases back, as the system ages. We show that the observed 
deviation from the Einstein relation is related to the relaxation times previously measured in dynamic 
light scattering experiments. 
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At thermodynamical equilibrium, the fluctuation- 
dissipation theorem (FDT) relates the response function 
of the system to its time autocorrelation function. This 
theorem may display several forms, depending on the 
considered observable (Nyquist formula, Einstein rela- 
tion), but involves a single thermodynamic parameter, 
which is the equilibrium temperature. However, FDT 
applies only to ergodic systems at equilibrium, and is 
not expected to hold for non-equilibrium systems, like 
glasses and gels, which exhibit relaxation times longer 
than or comparable to the observation timescale. 

Recently, many efforts have been devoted to apply sta- 
tistical physics concepts to such out-of-equilibrium sys- 
tems. In particular, fluctuation-dissipation relations have 
been extended with the help of a timescale dependent 
effective temperature, different from the bath tempera- 
ture, and that has been shown to display many of the 
properties of a thermodynamic temperature Q, Q ■ This 
deviation from equilibrium FDT has been observed in 

many numerical simulations [HHSEE0BIS[l3l- In 

these models and simulations, violations are expected to 
occur when the characteristic observation time 1/oj is of 
the same order or greater than the aging time t w . On 
the other hand, there are only few experimental studies 
of FDT violations in aging materials [3 13 13 13 
In all experimental systems, violations are observed when 
\/u> is smaller than t w . This suggests that the aging time 
t w is not the only relevant parameter to describe aging, 
as usually considered in numerical simulations, and that 
microscopic processes, charaterized by the distribution of 
relaxation times in the system, must be considered. 

In this letter we measure the evolution with aging time 
of an effective temperature T e g in a colloidal glass of 
Laponite. This is achieved by simultaneously measur- 
ing the mobility - using an optical tweezer - and the 
position fluctuations of a micrometric bead embedded in 
the glass. We find that, starting from the bath temper- 
ature, T e ff increases with t w , up to 1.8 times the bath 



temperature, and decreases back upon further aging. We 
propose a consistent interpretation of this unusual behav- 
ior by considering the evolution of the relaxation times 
distribution in the glass. 

Let us consider a diffusing particle of mass m evolv- 
ing in a stationary medium. Its motion is described 
by a generalized Langevin equation m dv/dt — F(t) — 
m "f(t — t') v(t') dt' , in which v(t) — dx/dt is the 
particle velocity, F(t) the random force acting on the 
particle and 7(f) a delayed friction kernel that takes into 
account the viscoelastic properties of the medium. If the 
surrounding stationary medium is in thermal equilibrium 
at temperature T, one can derive a generalized Einstein 
relation, which is a specific form of FDT, s 2 (Ax 2 (s)} = 
2kTfi(s), in which (Ax 2 (s)} is the mean-square displace- 
ment Laplace transform and fl(s) — l/m A f(s) the mobil- 
ity Laplace transform, if inertia is neglected. 

The general situation of a particle diffusing in an out- 
of-equilibrium environment is much more difficult to de- 
scribe. This point has been developed in details in [lfl | 
using a generalized Langevin formalism. The FDT is ex- 
tended with the help of a frequency-dependent effective 
temperature T e ff(w), parametrized by the age t w of the 
system, and which satisfies 0] : 



where 



T eS (uj)$ien(uj) = SJe[/x(u;)e(u;)] 



(Ax 2 {s)) = 2fc9(s) 



(1) 



(2) 



The two equations (1) and (2) define an out-of- 
equilibrium Einstein relation, for a given aging time t w . 
The parameter O(s) is related to Q(uj) — 0(s = — iuj). 

The principle of our experiment to measure the effec- 
tive temperature in an out-of-equilibrium system is the 
following : micrometric probes are immersed in a col- 
loidal glass to allow the measurements, at the same aging 
time t w , of both their fluctuating position and mobility. 
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T e g is then obtained from equations (1) and (2). The col- 
loidal glass consists in an aqueous suspension of Laponite 
RD, a hectorite synthetic clay provided by Laporte In- 
dustry. The preparation procedure of the glass has been 
addressed in details in [lq. These aqueous suspensions 
form glasses for low volume fraction in particles |l9( . Be- 
ing in a "liquid" state right after preparation, the sus- 
pension becomes more and more viscoelastic with time. 
Since the physical properties of the suspension depend 
on the time t w elapsed since preparation, the sample is 
said to age. Aging can be seen through the evolution 
of both the viscoelastic properties and of the colloidal 
disks diffusion [Isl Eoj]. This glass, obtained at the am- 
bient temperature, is optically transparent. Moreover, it 
presents other advantages that make possible to measure 
an effective temperature with tracer beads. First, the 
preparation procedure allows to obtain a reproducible 
initial state, leading to an accurate determination of the 
origin t w = 0. Second, the Laponite suspensions age on 
timescales that depend on the particles concentration. 
We are thus able to choose the aging timescales of the 
glass by adjusting this concentration. With a volume 
fraction of 2.3%wt, the glass evolves slowly enough to 
allow two successive measurements - fluctuation and dis- 
sipation - without significative aging of the sample. The 
two measurements are thus performed at the same t w . 

The experiments were carried out in a square chamber 
- 20 x 20 mm 2 - made of a microscope plate and a cover- 
slip separated by a thin spacer (0.1 mm thickness). The 
beads are suspended in the glass right after its prepa- 
ration. The chamber is then filled with the suspension, 
sealed with vacuum grease and mounted on a piezoelec- 
tric stage on the plate of an optical microscope. The 
probes are latex and silica beads, in very low concentra- 
tion (respectively 10 _4 % and 4.10 _4 % in volume). Latex 
beads (1.0±0.1/xm in diameter, Polysciences, Inc.), were 
preferentially used for fluctuation measurements : since 
they do not deposit during the experiment, their ran- 
dom motion is not perturbated by the chamber walls. 
Silica beads (2.1 ± 0.1/im in diameter, Bangs Lab Inc.) 
were used for dissipation measurements, because they are 
more efficiently trapped by the optical tweezers. The di- 
ameters of the two kinds of probes are close to each other, 
thus the comparison between the results of the fluctua- 
tion and dissipation measurements, once rescaled to the 
same diameter, is meaningful. 

Let us first focus on the two-dimensional brownian mo- 
tion of a tracer bead immersed in the glass. At a given 
aging time t w , we record the fluctuating motion of a 1 //m 
latex bead during 8 s, with a fast CCD camera sampling 
at 125 Hz (Kodak, PS-220). A digital image analysis al- 
lows to track the bead positions x(t) and y(t) close to 
the focus plane of the microscope objective. For each 
bead, we calculate the time-averaged mean-square dis- 
placement (Ar 2 {t)) t , = ((x{f +t)- x(t')) 2 + (y{f +t)- 
y(t')) 2 )t' — 2(Ax 2 (t)) t ' . To preserve a good statistics, 



we keep the data of (Ar 2 (t)) t > in the range 0.008 < t < 1 
s. The glass remains in a quasi-stationary state dur- 
ing the recording, which takes a short time compared to 
the aging timescale. The quantity (Ar 2 (t)) t i can thus 
be identified to the ensemble-averaged mean-square dis- 
placement {Ar 2 (t))E- Fig. 1 shows the mean-square 
displacement of latex beads immersed in the colloidal 
glass, as a function of time t, for various aging times 
t w . The mean-square displacement is well described by a 
power-law behavior (Ar 2 (t))E — Dt a over the full range 
0.008 < t < 1 s. Upon increasing on t w , the exponent 
decreases from 1.05 ± 0.05 at the earliest aging times 
to 0.25 ± 0.05 at long aging times. This indicates a 
nearly diffusive behavior of the tracer bead at short ag- 
ing times, that becomes sub-diffusive as the glass ages. 
We describe now the measurement of the mobility /u(o;), 
at a given frequency to for various aging times. We re- 
call that this measurement is performed right after the 
fluctuation motion recording, at the same aging time t w . 
As the Laponite suspension is a viscoelastic fluid, the 
bead mobility ^(w) = \iJ,(u))\e' l ^ ujS> is a complex num- 
ber. We thus need to measure the phase and modulus 
of the tracer mobility. We use an optical tweezer to trap 
a 2.1 fim silica bead immersed in the glass. Trapping 
is achieved by focusing a powerful infrared laser beam 
(Nd YAG, Spectra-Physics, P max = 600 mW) through 
a microscope objective of large numerical aperture [2 If . 
The trapping force F on a small dielectric object like a 
silica bead is proportional to the intensity gradient in the 
focusing region. To first order, one can write F = —kx, 
where x represents the distance of the trapped object 
from the trap center. The stiffness k is known from an 
independent calibration. Once the bead is trapped, we 
make the experimental chamber oscillate by monitoring 
the displacement x p expiiut) of a piezoelectric stage. As 
a consequence, the viscoelastic glass exerts a sinusodal 
force F' exp(iu>t) on the bead. We record with the fast 
camera the bead movement, and measure by conventional 
image analysis its displacement x exp (iu)t) from the trap 
center. At a given frequency uj, the force F'(oj) is given 




t(s) 



FIG. 1: Mean-square displacement of a 1 /im latex bead 
immersed in the glass, as a function of time. The curves cor- 
respond to different t w = 8, 125, 155, 178, 195 and 220 minutes 
from top to bottom. The fluctuating motion is purely diffu- 
sive at short t m and becomes sub-diffusive as the glass ages. 
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by F'(lu) — v(oj)/h(u>), where v(u>) = iu>(x p — x) is the 
relative glass / bead velocity, and fi(ui) the Fourier trans- 
form of the bead mobility. In our range of experimental 
frequencies (0.5 < / < 10 Hz), the bead inertia is negligi- 
ble, so that we can simply use the relation F + F' = to 
calculate |/x(w)| and <j>{u>). Notice that the piezoelectric 
stage and the camera are triggered by two synchronized 
signals, numerically generated by a PC computer, so that 
the phase shift between the force and the bead movement 
can be accurately measured. Fig. 2 (a,b) show the com- 
plex mobility of the tracer bead as a function of the aging 
time t w for various frequencies of the applied force. 

Two more steps are necessary to calculate the ef- 
fective temperature from both fluctuation and mobility 
data. First, the mean-square displacement is numerically 
Laplace transformed to the frequency domain 0.15 — 20 
Hz. In this range, we find that the Laplace transform 
is well adjusted by a power-law (Af 2 (s)) = as~ b . As 
expected from the observation (Ar 2 (t))E = Dt a , the 
relation b = a + 1 is accurately verified. Second, we 
analyze the frequency dependence of the mobility. In 
the experimental frequency range, is well fitted 

by a power law |^(w)| = fio^ , a s shown in Fig. 2c. 
The exponent f3 increases with t w from zero at low ag- 
ing times to 0.75 ± 0.05 at the end of experiment. To 
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FIG. 2: Modulus (a) and phase (b) of the complex mo- 
bility fi(u>) = \fi(u})\e"^^ of the tracer bead as a function 
of t w , for various frequencies of the applied force. The full 
circles in (b) correspond to (j> = — /3(t w )n/2; (c) : mobil- 
ity modulus versus ui for various t w ; from top to bottom, 
t w — 120, 145, 165, 183, 205 min. At low t w , is indepen- 

dent of lo. Upon increasing t w , the modulus is well fitted by 
a power-law = Hou 13 with f3 only depending on t w - 



a first-order approximation, we consider that the phase 
(f> is independent of the frequency. Within a good ap- 
proximation, its dependence on t w can be related to (3 by 
cf> = —f3(t w )ir/2 as shown in Fig. 2b. From the analytical 
form n(u)) = exp (— i/?7r/2), the Laplace transform 
/i(s) = fi(u> = is) = /.los^ is derived by analytical continu- 
ation. Using equation (2), we calculate the function O(s), 
parametrized by t w , and derive <d(u>) = 9(s = — iui). Fi- 
nally, the effective temperature, at a given t w , is : 



kT eS (u) 



a cos((b-2)n/2) 
4// cos(Pn/2) 



,2-6-/3 



(3) 



The dependence of T e g on t w is shown in Fig. 3, at 
a frequency / = 1 Hz. The results have been averaged 
over three realizations. At the earliest t w , the effective 
temperature is close to the bath temperature Th a th = 300 
K. Upon increase on t w , T c g increases up to 1.8 times 
the bath temperature and then decreases towards Tt, a th 
upon further increase on t w . This is the first time that 
such a behavior - increase of T e g followed by a decrease 
- is observed in a glassy system. A realistic explanation 
of this non-monotonic behavior is proposed in the next 
paragraph. Moreover, our results are in contradiction 
with electrical measurements of FDT performed in the 
same system 0] , where T e g is found to decrease with t w 
and w and is larger by about one order of magnitude. A 
possible origin of this discrepancy may be the choice of 
different observables in the two experiments [Toj| . 

An interpretation of the dependence of T e g on t w is 
provided by Dynamic Light Scattering (DLS) and Diffu- 
sive Wave Spectroscopy (DWS) experiments, previously 
performed in colloidal glasses of Laponite |2tJ, |22| . The 
resulting distribution function of relaxation times P(t) is 
schemed in Fig. 4. Upon increasing t w , part of the modes 
distribution function shift to larger times, while the mode 
at t ~ 0.1 ms remains unchanged. When probing our 
experimental system at a typical frequency / = 1 Hz, 
three situations occur. At the earliest t w , the fast modes 
t << 1// allow the system to thermalize with the bath 
(Fig. 4a). The measured effective temperature is the 
bath temperature. Upon increase on t w , modes r ~ 1// 
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FIG. 3: Effective temperature of the colloidal glass versus tw, 
measured at a frequency / = 1 Hz. Upon increase on t w , T c g 
increases up to 1.8 times the bath temperature and decreases 
back to Tbath- The dashed line is a guide for the eyes. 
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appear in the system as seen in Fig. 4b. On this observa- 
tion timescale 1//, the system is out-of-equilibrium : the 
measured temperature becomes different from the bath 
temperature and is timescale dependent. Deviation from 
the equilibrium Einstein relation is thus observed. Fi- 
nally, for very long t w , the fast modes allow the system 
to thermalize with the bath while the slow ones t>1// 
do not play any role at the experimental timescale. T c s 
is then expected to reduce back to Tbath- These different 
situations are clearly identified in Fig. 3. However the 
experimental set-up does not allow to observe a further 
decrease of T e g at long t w . Indeed, beyond t w = 205 min- 
utes, the mobility modulus becomes smaller than 10~ 2 m 
s _1 pN _1 . In this range, the optical tweezer is not pow- 
erful enough to induce a detectable motion of the bead. 

In these experiments, deviations are observed when 
lot ~ 1. This suggests that, besides the waiting time 
t w , the distribution of relaxation times must be included 
in models to achieve an accurate description of aging. 

Let us now comment the dependence of T e g on the 
frequency. Since we experimentally find that 2 — b — (3 > 
at all t w , T ff increases with to as w 2 ~ & ~^, according 
to equation (3). This behavior is consistent with our 
interpretation at least in the first situation. However, 
in our picture, one expects T c g to decrease with u> at 
long t w . Actually, one can see from Fig. 2b that <f> can 
no longer be accurately considered as independent of w 
at long t w . Thus, the analytical interpolation fJ,(u>) — 
[i^io^ exp — iPtt/2 probably breaks down. This difficulty 
will be soon overcomed by a direct measurement of /t(s) 
from a creep experiment. 

This work provides a test of the nonequilibrium Ein- 
stein relation in a colloidal glass, using diffusion and mo- 
bility measurements on a micrometric probe. We observe 
for the first time in a glassy system a non-monotonic be- 
havior of T e ff with t w . It seems likely that this behavior 
is directly related to the evolution of the relaxation times 
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FIG. 4: Scheme of the distribution function of relaxation 
times P(t) in the glass of Laponite (typically 2.5% wt) at dif- 
ferent t w . Upon increasing t w , part of the modes distribution 
shifts towards larger times, while the mode at r ~ 0.1 ms 
remains unchanged. The arrow represents the measurements 
timescale 1//. 



distribution. Notice that deviations from the Einstein re- 
lation, found in supercooled liquids even at equilibrium, 
were explained by spatial heterogeneities |23j. Future 
works will have to analyze how such heterogeneities in- 
terplay with the concept of T e ff - which is an ensemble- 
averaged quantity. Another open question is whether this 
effective temperature has a real thermodynamics mean- 
ing plllflj. The answer will come from experimental tests 
of FDT involving other physical observables. 
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